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Abstract 

We present a new short-time approximation scheme for evaluation of decoherence. 
At low temperatures, the approximation is argued to apply at intermediate times as 
well. It then provides a tractable approach complementary to Markovian-type approx- 
imations, and appropriate for evaluation of deviations from pure states in quantum 
computing models. 
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1. Introduction 

Consider a microscopic quantum system with the Hamiltonian H$- We will refer 
to the quantum-computing single quantum bit (qubit) or multi-qubit paradigm to help 
define the questions and set up the challenges, in describing how the system, S, interacts 
with the surrounding macroscopic world. However, in principle S can be any quantum 
system. 

Interactions with the surroundings can be quite different depending on the setting. 
For example, in quantum measurement, which for orthodox quantum theory is not fully 
understood, the wavefunction of the system is probed, so part of the process would 
involve a strong interaction with the measuring device, such that the system's own 
Hamiltonian plays no role in the process. However, in most applications, the external 
interactions are actually quite weak. Furthermore, the aim is to minimize their effect, 
especially in quantum computing. 

Traditionally, interactions with the surrounding world have been modeled by the 
modes of a bath, B, with each mode described by its Hamiltonian Mjf, so that the bath 
of modes is represented by 

H b = J2 m k- (1-1) 

K 

The interaction, I, of the bath modes with the system S, will be modeled by 

H i = A s P b =AsJ2 J k > ( L2 ) 

K 

where As is some Hermitean operator of S, coupled to the operator Pb of the bath. 

The bath, or "heat bath" , can be a collection of modes, such as photons, phonons, 
spins, excitons, etc. For a bosonic bath of oscillators, [1-6], which we use for derivation 
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of specific results, we take 



M K = u K a ] K aK , (1-3) 

Jk = g*K a K + gKa ] K . (1.4) 

Here we have assumed that the energy of the ground state is shifted to zero for each 
oscillator, and we work in units such that h = 1. 

The total Hamiltonian of the system and bath is 

H = Hs + Hb + Hj . (1.5) 

More generally, the interaction, (1.2), can involve several system operators, each cou- 
pling differently to the bath modes, or even to different baths. The bath modes, in turn, 
can be coupled to specified external objects, such as impurities. 

Let p(t) represent the reduced density matrix of the system at time t > 0, after 
the bath modes have been traced over. For large times, the effect of the environment 
on a quantum system that is not otherwise externally controlled, is expected to be 
thermalization: the density matrix should approach 

/. \ exp(-f3H s ) 

" {t ^ °°> = 5JSFw| • (L6) 

where (3 = 1/kT. At all times, we can consider the degree to which the system has 
departed from coherent pure-quantum-state evolution. This departure is due to the 
interactions and entanglement with the bath. We also expect that the temperature, T, 
and other external parameters that might be needed to characterize the system's density 
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matrix, are determined by the properties of the bath, which in turn might interact with 
the rest of the universe. 

Let us introduce the eigenstates of Hs, 

H s \n} = E n \n), (1.7) 

and have AE denote the characteristic energy gap values of S. We also consider the 
matrix elements of p(t), 

Pmn(t) = (m\p(t)\n) . (1.8) 

For large times, we expect the diagonal elements p nn to approach values proportional 
to e~ liEn , while the off-diagonal elements, p m / n , to vanish. These properties can be 
referred to as thermalization and decoherence in the energy basis, though "thermaliza- 
tion" in the strong sense of (1.6) implies decoherence. 

To establish these thermalization and decoherence properties, several assumptions 
are made regarding the system and bath dynamics [1-11]. At time t = 0, it is usually 
assumed that the bath modes, K, are thermalized, i.e., have density matrices 

9 K = e ~^ MK /Tv k (e-? MK ) . (1.9) 

The density matrix R of the system plus bath at time t = is assumed to be the direct 
product 

R(o) = p(o)Y[e K , (i.io) 

K 

and the system and bath modes are not entangled with each other. 
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Now, a series of assumptions are made, e.g., the Markovian and secular approxima- 
tions. The most important is the Markovian approximation, which, even though it can 
be stated and introduced in various ways, essentially assumes that the density matrices 
of the bath modes are reset externally to the thermal ones, on time scales shorter than 
any dynamical times of the system interacting with the bath, and the product form of 
the full density matrix is maintained. This is a natural assumption, because each bath 
mode is coupled only weakly to the system, whereas it is "monitored" by the rest of 
the universe and kept at temperature T. In its straightforward version, this amounts to 
using (1.10) for times t > 0. Ultimately, such approaches aim at master equations for 
the evolution of p m n(t) at large times, consistent with the Golden Rule and with the 
expected thermalization and decoherence properties. 

In variants of these formalisms, several time scales are identified. One is the inverse 
of the upper cutoff, Debye frequency of the bath modes, 1/ujd- Another is the thermal 
time h/kT = (3 (in units of fi = 1). The system S has its own characteristic time, 
1/AE, as well as the system-bath dynamical times of thermalization and decoherence, 
etc., Ti )2 ,..., corresponding to the "intrinsic" NMR/ESR times Ti, T 2 , etc. Heuristically, 
bath modes of frequencies u> comparable to AE are needed to drive thermalization and 
decoherence. Initial decoherence can be also mediated by the modes near to = 0. At 
low temperatures, we can assume that 1/ud < 1/AE < (3. 

There is evidence [7,11,12] that at low temperatures, the Markovian- type and other 
approximations used in the derivation of equations for thermalization and decoherence, 
are only valid for times larger than the thermal time scale (3. For quantum comput- 
ing applications, in solid-state semiconductor-heterostructure architectures [13-19], we 
expect temperatures of several tens of //K. The thermal time scale then becomes dan- 
gerously close to the external single-qubit control, Rabi-flip time even for slower qubits, 
those based on nuclear spins. We emphasize that not all the approximation schemes 
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have this problem [11]. 

In Section 2, we offer additional comments on decoherence and quantum comput- 
ing. Then, in Section 3, we develop a short-time-decoherence approximation. In a 
discussion at the end of Section 3, we offer arguments that, at low temperatures, our 
approximation is actually valid for intermediate times, larger than 1/ud, hopefully up 
to times comparable or larger than 1/AE. Specific results for the bosonic heat bath are 
presented in Section 4. Section 5 comments on the case of adiabatic decoherence, when 
the short-time approximation becomes exact. 
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2. Decoherence and quantum computing 

Quantum computing architectures usually emphasize systems, both the qubits and 
the modes that couple them (and at the same time act as a bath mediating unwanted 
coupling to the rest of the universe), that have large spectral gaps. It is believed that, 
especially at low temperatures, spectral gaps slow down relaxation processes. Therefore, 
quantum computing architectures usually assume [13-19] qubits in quantum dots, or in 
atoms, or subject to large magnetic fields, and coupled by highly nondissipative quantum 
media [14,19]. 

The spectral gaps are expected to slow down exponentially, by the Boltzmann 
factor, the processes of thermalization, involving energy exchange. Off-shell virtual 
exchanges, will be also slowed down, but less profoundly. The latter processes contribute 
to decoherence. Therefore, at low temperatures, we might expect separation of time 
scales of the initial decoherence vs. later-stage thermalization and further decoherence. 
The latter two processes are described by the traditional NMR/ESR intrinsic T\ and 
T 2 , respectively. 

Since only thermalization is clearly associated with the energy eigenbasis, one can 
also ask whether the energy basis is the appropriate one to describe decoherence for 
short and intermediate times, before the thermalizing processes, that also further drive 
decoherence, take over. The issue of the appropriate basis for studying decoherence, has 
also come up in models of quantum measurement. It has been argued [20-24] that the 
eigenbasis of the interaction operator, As, may be more appropriate for intermediate 
times than the energy eigenbasis. 

Yet another aspect of decoherence in quantum computing, involves the observation 
that we really want to retain a pure state in the quantum computation process [25-30] . 
Decay of off-diagonal matrix elements, in whatever basis, might not be the best measure 
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of deviations from the pure-state density matrix, where by pure states we mean those 
with density matrices that are projection operators For instance, the deviation 

of Tr s [p 2 (t)] from 1, may be more appropriate, and is easer to calculate than other 
measures, specifically those motivated by the "entropic" expressions proportional to 
Tr s [pln(p)] . Therefore, it is desirable to have basis-independent expressions for the 
reduced density operator p(t). 

Recently, several groups have reported [12,19,24,31-41] results for spin decoherence 
in solid state systems appropriate for quantum computing architectures. Some of these 
works have not invoked the full battery of the traditional approximations, Markovian 
and secular, etc., or have utilized the spectral gap of the bath modes, to achieve better 
reliability of the short-time results. In [41], interaction of the spin-exciton bath modes 
with impurities was accounted for, as the main mechanism of decoherence. In the 
present work, we limit ourselves to the bath modes only interacting with the system. 
Experimental efforts are picking up momentum, with the first limited results available 
[42,43] by traditional NMR/ESR techniques, with the quantum-computing emphasis. 

An approach, termed adiabatic decoherence, have been developed by us [24], ex- 
panding the earlier works [12,31-33], with the goal of avoiding the ambiguity of the basis 
selection and achieving exact solvability. The price paid was the assumption that Hs 
is conserved (a particular version of the quantum nondemolition processes), which is 
equivalent to requiring that 

[H s , H] = [H s , As] = (adiabatic case) . (2.1) 

This makes the eigenbasis of Hs and As the same, but precludes energy relaxation, 
thus artificially leaving only energy-conserving relaxation pathways that contribute to 
decoherence. We will comment on the results of this approach in Section 5. 
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Most of the results referred to earlier, have involved approximations of one sort or 
another. The most popular and widely used approximation has been the second-order 
perturbative expansion in the interaction strength, if/, though some nonperturbative 
results have also been reported. In Section 3, we describe a novel approximation scheme 
[44] that is valid for short times. It has several advantages, such as becoming exact in 
the adiabatic case, allowing derivation of several explicit results, and, at least in princi- 
ple, permitting derivation of higher-order approximations. Certain models of quantum 
measurement evaluate decoherence by effectively setting Hs = 0. Our approximation 
then becomes exact, and our results are consistent with these studies [45,46]. 

Our formulation in Section 3, will be quite general, and we will not use the specific 
bath or thermalization assumptions. However, we do utilize the factorization property 
(1.10) at time t = 0. Thus, we do have to assume that, at least initially, the system and 
the bath modes are not entangled. In fact, the present formulation also relies on that 
the Hamiltonians at hand are all time-independent. Therefore, we have excluded the 
possibility of controlled dynamics, in the quantum computing sense, when gate functions 
are accomplished by external couplings to individual qubits and by external control of 
their pairwise interactions. Our formulation, therefore, applies to "idling" qubits or 
systems of (possibly interacting) qubits. It is reasonable to assume that a lower limit 
on decoherence rate can be evaluated in such an idling state, even though for quantum 
error correction, qubits otherwise idling, might be frequently probed (measured) and 
entangled with ancillary qubits [25-30]. 

The t = factorization assumption (1.10), shared by all the recent spin-decoherence 
studies, then represents the expectation that external control by short-duration but large 
externally applied potentials, measurement, etc., will "reset" the qubits, disentangling 
them from the environment modes to which the affected qubits are only weakly coupled. 
Thus, we assert that it is the qubit system that gets approximately reset and disentan- 
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gled from the bath towards time t = 0, instead of the bath being thermalized by the 
rest of the universe, as assumed in Markovian approximation schemes. 
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3. Short-time decoherence 

In addition to the energy basis, (1.7), we also define the eigenstates of the interaction 
operator As, by 



As|7) = A 7 |7> , 



(3.1) 



where the Greek index labels the eigenstates of As, with eigenvalues A 7 , while the 
Roman indices will be used for the energy basis, and, when capitalized, for the bath 
modes, (1.2)-(1.4). 

The time dependence of the density matrix R(t) of the system and bath, is formally 
given by 



We will utilize the following approximate relation for the exponential factors, as our 
short-time approximation, 



This relation has the following appealing properties. It becomes exact for the adiabatic 
case, (2.1). Furthermore, if we use the right-hand side and its inverse to replace e ±lHt , 
then we are imposing three time-evolution-type transformations on R(0). Therefore, 
the approximate expression for R(t) will have all the desired properties of a density 
operator. Finally, extensions to higher-order approximations in powers of t are possible, 
by using relations derived in [47], where various expressions valid to 0(t 4 ) and 0(t 5 ) 
were considered. 




(3.2) 



^(Hs+HB+H^t+Oit 3 ) _ e iH s t/2 e i(HB+H!)t e iH s t/2 



(3.3) 



Our goal is to evaluate the resulting approximation to the matrix element, 



Pmn (t) = ^ B ^ m \ e -^st/2 e ^(H B+Hl )t e ^H s t/2 R ^^Hst/2 ^(H B+Hl )t ^H s t/2^ n) 

(3.4) 

First, we apply the operators H$ in the outer exponentials, acting to the left on (m| , 
and to the right on \n), replacing H$ by, respectively, E m and E n . We then note that 
the second exponential operator in (3.4) contains As, see (1.2). Therefore, we insert the 
decomposition of the unit operator in the system space, in terms of the eigenbasis of As, 
before the second exponential, and one in terms of the eigenbasis of Hs after it. This 
allows us to apply As in the second exponential and also Hs in the third exponential. 
The same substitution is carried out on the other side of -R(O), with the result 



Pmn{t) = ^ Tr B 
7 p q 8 



e -^V2 (m | 7)(7 | p)e -i(H B +A T P B )t e -^t/2 p ^ (0) 



K 



iE q t/2 i(H B +\ s P B )t 



(q\S){8\n)e iE ^ 2 ] . 



(3.5) 



The next step is to collect all the terms, and also identify that the trace over the 
bath can be now carried out for each mode separately. We use (1.1)-(1.2) to write 



.(«) = £ {< 



i(Eq+E n — E p — E m )t/2 



(m| 7 ><7b)p M (0) <#><%> 



x JJ rj^ K ^ e -i(M K +\, / J K )t k e i(M K +X s J K )t | _ 
K 



(3.6) 



While this expression looks formidable, it actually allows rather straightforward calcula- 
tions in some cases. Specifically, the simplest quantum-computing applications involve 
two-state systems. Then the sums in (3.6) are over two terms each. The calculations 



involving the overlap Dirac brackets between the eigenstates of H$ (labeled by m, n, 
p and q) and of Ag (labeled by 7 and S), as well as the energy-basis matrix elements 
of p(0), cf. (1.8), involve at most diagonalization of two- by- two Hermitean matrices. 
Of course, the approximation (3.6) can be used for evaluation of short-time density 
matrices for systems more general than two-state. 

The challenging part of the calculation involves the trace over each mode of the 
bath. Since these modes have identical structure, e.g., (1.3)-(1.4) for the bosonic bath 
case, but with X-dependent coupling constants, the calculation needs only be done 
once, in the space of one mode. Furthermore, results for the bath models ordinarily 
used, such as the bosonic and spin baths, are either already available in the literature 
or can be calculated without much difficulty. For the thermalized initial bath-mode 
density matrix 6*^, we give the exact bosonic- model expression in the next section. 

In the remainder of this section, we first further analyze the trace over one bath 
mode entering (3.6). We then comment on the limits of validity of the present approx- 
imation. 

In an obvious shorthand notation, we write the single-mode trace in (3.6) as 



^ e -i(M+~/J)t q e i(M+5J)t 



Tr Q e i i M + 5J ) t e -i(M+j.J)t 



(3.7) 



Now, to the same order of approximation as used in (3.3), we can write 



i(M+5J)t+0(t 3 ) _ iMt/2 iSJt iMt/2 



(3.8) 



The resulting approximation for the trace (3.7) reads 




(3.9) 



which illustrates that, within this approximation, the product of traces in (3.6) is a 
function of the difference A 7 — X§- In fact, this product is exactly 1 for A 7 = \$ and, in 
most applications, the following form is likely to emerge, 



n^H ] 



- const (X^-X s ) 2 t 2 +0(t 3 ) 



(3.10) 



K 



though we caution the reader that (3.10) is somewhat speculative and suggested by the 
exact result for the bosonic heat bath, reported in the next section. 

Finally, we point out that in most cases of interest, the initial single-mode density 
matrix 9 will commute with the bath- mode energy operator M. In fact, the thermalized 
is a function of M. Therefore, (3.9) can be further simplified to 



However, let us emphasize that the approximate relations (3.9)-(3.11) are likely of 
value only as far as they help to derive basis-independent (operator) approximations to 
p(t), by a technique illustrated in the next section. Indeed, for most bath models it is 
advisable to calculate the single-mode trace exactly first, according to (3.6), and then 
attempt various approximations. 

The latter statement reflects our expectation that the approximation developed 
here is valid, for low temperatures, not only for short times, defined by t < 1/u>d, 
but also for intermediate times, exceeding 1/ujd- This is suggested by the result of 
an illustrative calculation in the next section, but mainly by the fact that (3.11) only 
includes the bath-mode energy scales via 6>, and, therefore, at low temperatures, is 
dominated by the lowest bath-mode excitations, and is not sensitive to frequencies 
of order ud- Thus, we expect our approximation to be applicable complementary to 



Tr Oe^ 5 -^ 



(3.11) 



the Markovian-type approximations and definitely break down in the regime of fully 
developed thermalization, for t > 0((3). Additional supporting observations are offered 
in Section 5, when we consider the adiabatic case (2.1). 
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4. The bosonic heat bath 

In this section, we consider the bosonic heat bath [6], see (1.3)-(1.4), in the initially 
thermalized state, 



9 K = e- /3MK /Tr K (e- f3MK ) = (i- e > K )e-^ a ^ x . (4.1) 

The product of the single- mode traces in (3.6), is then available in the literature 
[12,24,31], 



Pmn{t)= {^ Eq+En - Ep - Em)t/2 {m\l){l\p){q\5){5\n)p P M 



jp qS 



K 



x exp ( - J2 ^jf~ [ 2 ( A 7 - A *) 2 sin2 coth ^ + i W - X ") ( sin 00 Kt ~ ^ } ) } 



(4.2) 



The last term in the exponent, linear in t, is usually viewed as "renormalization" of the 
system energy levels due to its interaction with the bath modes. It can be removed by 
adding the term, 



H r = A 2 s Y,\9k\ 2 /uk, (4-3) 

K 

to the total Hamiltonian. However, the usefulness of this identification for short times 
is not clear, and we will not use it. One can check that, unmodified, (4.2) is consistent 
with the expectation (3.10). 

Let us now define two no n- negative real spectral sums, B(t) and C(t), over the 
bath modes, 
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B\t) 



8 



E 



K 




coth 



/3u K 
2 



(4.4) 



C(t) 



E 



K 




(4.5) 



When converted to integrals over the bath mode frequencies, with the cutoff at lod, 
these sums have been discussed extensively in the literature [6,12,31], for several choices 
of the bath mode density of states and coupling strength g as functions of the mode 
frequency. 

The final expression is, 



When the spectral functions are expanded in powers of t, this result confirms all the 
conclusions and conjectures discussed in Section 3, in connection with relations (3.9)- 



Let us now turn to the derivation of the basis-independent representation for p(t), 
by utilizing the integral identity 



Exponential factors in (4.6) can then be reproduced by applying operators on the wave- 
functions entering the overlap Dirac brackets, with the result 





(4.6) 



(3.11). 




(4.7) 
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yfifp(t) = Jdye~ y2 e~ lHst/2 e *[wS(t)A s +C(t)A|] e ~iH s t/2 p ^ e iH s t/2 e -i[yB(t)A s +C(t)A 2 s ] e iH s t/2 

(4.8) 

Within the 0(t 2 ) approximation (3.3), given that B and C are of order linear or 
higher in t, we can combine the exponential operators to get an alternative approxima- 
tion, 




though (4.6) and (4.8) are in fact easier to handle in actual calculations. 

As an application, let us consider the case of H$ proportional to the Pauli matrix 
cr 2 , e.g., a spin-1/2 particle in magnetic field, and As = cr x , with the proportionality 
constant in the latter relation absorbed in the definition of the coupling constants gx 
in (1.4). Let us study the deviation of the state of a spin-1/2 qubit, initially in the 
energy eigenstate | f) or | |), from pure state, by calculating Tr s [p 2 (t)] according to 
(4.8). We note that for a two-by-two density matrix, this trace can vary from 1 for pure 
quantum states to the lowest value of 1/2 for maximally mixed states. 

A straightforward calculation with p(0) = |T)(tl or li)(i|> yields 



1 + e 



-2B 2 {t) 



(4.10) 



As the time increases, the function B 2 (t) grows monotonically from zero [6,12,24,31]. 
Specifically, for Ohmic dissipation, B 2 (t) increases quadratically for short times t < 
0(1/u>d), then logarithmically for 0(1/ud) < t < 0(h/kT) : and linearly for t > 
0(h/kT). For other bath models, it need not diverge to infinity at large times. 



Both approximations, (4.8)-(4.9), make the deviation from a pure state p(0) = 
|V'o)(' i /'o| apparent: p(t > 0) is obviously a mixture (integral over y) of pure-state pro- 
jectors \ip(y,t))(ip(y,t)\, where, for instance for (4.9), 

t) = e -i[tH S -yB(t)A s -C(t)A%] ^ ^ (4 _ n) 

with a somewhat different expression for (4.8). 
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5. The adiabatic case 

Relation (2.1) corresponds to the system's energy conservation. Therefore, energy 
flow in and out of the system is not possible, and normal thermalization mechanisms 
are blocked. This "adiabatic decoherence" limit thus corresponds to "pure dephasing" ; 
see [48]. 

The fact that our approximation becomes exact in this case, provides support to 
the expectation that, at low temperatures, it is generally valid beyond the cutoff time 
scale 1/u>d, providing a reasonable evaluation of decoherence and deviation from a pure 
state, as exemplified by the calculation yielding (4.10), in Section 4. 

With (2.1), we can select a common eigenbasis for H$ and Ag. Then the distinction 
between the lower-case Roman and Greek indices in (3.6) becomes irrelevant, and the 
sums can all be evaluated to yield 



Pmn(t) = e l ^- Em)t p mn (0) ]\ J e -ii M K+\mJK)t q k e i{M K+ \ n J K )t\ _ (5 ^ 



K 



This expression was discussed in detail in our work on adiabatic decoherence [24] . Specif- 
ically, for the initially thermalized bosonic heat bath case, we have, for the absolute 
values of the density matrix elements, 



|pm»(f)| = \pmn(0)\e- B2 ^ X -- X ^^ . (5.2) 

The decay of the off-diagonal matrix elements thus depends of the properties of the 
spectral function B 2 (t) as the time increases. Such explicit results [12,24,31-33] illustrate 
that for irreversibile behavior, the number of bath modes must be infinite, with the 
spectral function evaluated in the continuum limit. 
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In summary, we have derived short-time approximations for the density matrix and 
its energy-basis matrix elements. Our expressions are quite easy to work with, because 
for few-qubit systems they only involve manipulation of finite-dimensional matrices, and 
they will be useful in estimating decoherence and deviation from pure states in quantum 
computing models, including results for low temperatures. 

This research was supported by the National Science Foundation, grants DMR- 
0121146 and ECS-0102500, and by the National Security Agency and Advanced Re- 
search and Development Activity under Army Research Office contract D A AD- 19-99- 
1-0342. 
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